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Abstract. We introduce a new geometric evolution equation for hy- 
persurfaces in asymptotically flat spacetime initial data sets, that unites 
the theory of marginally outer trapped surfaces (MOTS) with the study 
of inverse mean curvature flow in asymptotically flat Riemannian man- 
ifolds. A theory of weak solutions is developed using level-set methods 
and an appropriate variational principle. This new flow has a natural 
application as a variational- type approach to constructing MOTS, and 
this work also gives new insights into the theory of weak solutions of 
inverse mean curvature flow. 

In what follows we consider an initial data set (M n+1 , g, K) that arises 
as a spacelike hypersurface M n+1 in a Lorentzian spacetime (L n+2 , h), with 
induced metric g and second fundamental form tensor K. We further as- 
sume that the initial data set (M,g,K) is asymptotically flat, that is, there 
exists a compact set O C M such that M\Q consists of a finite number of 
components, each diffeomorphic to R n+1 \.B(0, 1) and such that under these 
diffeomorphisms, the metric tensor g and second fundamental form K satisfy 

. . C C 
\9ij - d ij\ ^ mTT^T' \9ij,k\ < rn^, 
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as |x| —> oo, where the derivatives are taken with respect to the Euclidean 
metric. 

Let n denote the future directed timelike unit normal vector field of M C 
L, and consider a 2-sided hypersurface S n C M n+1 with globally defined 
outer unit normal vector field v in M. The mean curvature vector of S 
inside the spacetime L is then given by 

# E := Hv- Pn, 

where H := div^fV) denotes the mean curvature of £ in M, and P := tr-^K 
is the trace of K over the tangent space of S. 

The new initial value problem is then defined as follows. Given a smooth 
hypersurface immersion Fq : S — > M, the evolution of So := Fq(E) by in- 
verse null mean curvature is the one-parameter family of smooth immersions 
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F : E x [0, T) -)• M satisfying 

(*) (f ( ^ = iffp (M) ' ^E,*>0, 

1 F(-,0) = F . 

The quantity H+P corresponds to the null expansion or null mean curvature 
9j, of E t := F(E, f) with respect to its future directed outward null vector 
field l + := v + ft, 

6+ := (Hx t ,l+) h = H + P, 
and we assume that {H + -P)| E > so that (*) is parabolic and the surface 
Tif expands under the flow. This flow is a generalisation of inverse mean 
curvature flow, which corresponds to the special time-symmetric case of (*) 
where K = 0. Analogous to inverse mean curvature flow, in general it is 
expected that the null mean curvature of solutions of (*) will tend to zero at 
some points, and that singularities will develop. Therefore, the main part of 
this work is devoted to developing a theory of weak solutions of the classical 
flow (*). 

The motivation for introducing this particular generalisation of inverse 
mean curvature flow follows from the study of black holes and mass/energy 
inequalities in general relativity. In particular, it is hoped that this new flow 
will help to gain insight into the long standing Penrose conjecture in general 
relativity, which generalises the Riemannian Penrose inequality, proven by 
Huisken and Ilmanen [10] using their theory of weak solutions to inverse 
mean curvature flow, (see [3] for an alternative proof by Bray, which applies 
to the case of multiple horizons). 

More specifically, this flow is motivated by the theory of marginally outer 
trapped surfaces in general relativity. Physically, the outward null mean 
curvature 9^ measures the divergence of the outward directed light rays 
emanating from E. If 9^ vanishes on all of E, then E is called a marginally 
outer trapped surface, or MOTS for short. MOTS play the role of apparent 
horizons, or quasi-local black hole boundaries in general relativity, and are 
particularly useful for numerically modelling the dynamics and evolution of 
black holes. 

From a mathematical point of view, MOTS are the Lorentzian analogue of 
minimal surfaces. However, since MOTS are not stationary solutions of an 
elliptic variational problem, the direct method of the calculus of variations 
is not a viable approach to the existence theory. One successful approach to 
proving existence of MOTS comes from studying the blow-up set of solutions 
of Jang's equation 

I vw» \ I v.v^ + \ _ 

V |V«f + V VVlVuf + l ') 

for the height function wofa hypersurface, which was an essential ingredient 
in the Schoen-Yau proof of the positive mass theorem [15]. In their analy- 
sis, Schoen and Yau showed that the boundary of the blow-up set of Jang's 
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equation consists of marginally trapped surfaces. Building upon this work, 
existence of MOTS in compact data sets with two boundary components, 
such that the inner boundary is (outer) trapped and the outer boundary 
is (outer) untrapped, has been proven by Andersson and Metzger [2], and 
independently by Eichmair [6] , using different approaches for the barrier ar- 
gument at the boundary. Similarly, we see below that Jang's equation also 
plays a key role in the existence theory for weak solutions of (*). 

To develop the weak formulation for the classical evolution (*), we use the 
level-set method and assume the evolving surfaces are given by the level-sets, 

(2) S t = d{x e M | u(x) < t}, 

of a scalar function u : M — )■ M. Then whenever u is smooth and Vu ^ 0, 
the surface flow equation (*) is equivalent to the following degenerate elliptic 
scalar PDE 

V«\ . / ,-„• V'uV^u 



<**> d ""{m) + { 3 "-i^r) K " = ^ ul 

However, since this is an expanding flow, the function u is monotone non- 
decreasing. Therefore it only makes sense to study (**) on initial data sets 
(M,g,K) satisfying tvj^K > outside So, so that the zero function is a 
subsolution barrier for the Dirichlet problem for (**). 

In order to solve (**), we employ the method of elliptic regularisation, 
and study solutions, u £ , of the following strictly elliptic equation 

Vu r \ ( Si VVV-V 



A notable feature of elliptic regularisation that is heavily exploited in this 
work is that the downward translating graph 

'u F t" 



(3) Sf:= graph ^ 

solves the classical evolution (*) in the product manifold (M x M.,g := 
g®dz 2 ), where we extend the given data K to be parallel in the z-direction. 
Furthermore, this elliptic regularisation problem sheds new light on the 
study of Jang's equation ([T]), since a rescaling of (**) e can be interpreted as 
([T|) with a gradient regularisation term. 

To define weak solutions to (**), we use a variational principle for the 
energy functional 

(4) J U A U (F) := \d*F n A\ - [ \Vu\ - (j* - ^V) K ij} 

J FnA 

defined for sets F of locally finite perimeter and any compact set A. Here 
d*F denotes the reduced boundary of F, and v represents the unit normal 
Vu/\Vu\ to the surfaces S t defined by ([2]). The special case K = cor- 
responds to the functional employed by Huisken and Ilmanen in |10| . and 
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weak solutions of (**) necessarily exhibit the same jumping phenomenon, 
characteristic of weak solutions of inverse mean curvature flow. However, 
since Vtt/|Vu| is undefined on plateaus of the locally Lipschitz function u, 
we must define an appropriate notion of normal vector in these jump re- 
gions. For this reason, a careful analysis of the jump region of the limit of 
the elliptic regularisation problem (*) £ is vital to determining the correct 
formulation of weak solutions to (*), and constitutes a significant part of 
this work. On the contrary, a complete analysis of jump regions of inverse 
mean curvature flow is not included in }10| . since it was not necessary for 
the proof of the Riemannian Penrose Inequality. 

The main result of this work is the following weak existence theorem. 

Theorem 1. Let (M,g,K) be a complete, connected, asymptotically flat 
initial data set and let Eq be any precompact, smooth open set in M . Assume 
that trMK > on M\Eq. Then there exists a weak solution of (**) in M 
with initial condition Eq. 

The variational principle defining weak solutions also leads to a geometric 
characterisation of the flow, and in particular, of jump regions. We show 
that the level sets Et are outward optimising (see (|55f) ) in the sense that 
they minimise "area plus bulk energy P" on the outside, along the family of 
surfaces. This one-sided variational principle can then be exploited, via the 
choice of an appropriate initial condition, to prove the following existence 
result for MOTS. 

Proposition 2. The weak solution of (**) with outer trapped initial con- 
dition Eq satisfying 6q Eq < will jump immediately to a smooth MOTS in 
M\E at time t = 0. 

Proposition [2] highlights the natural utility of this flow as a variational 
type approach to constructing MOTS in initial data sets containing an outer 
trapped surface Eq, where trpfK > on M\Eq. 

We remark that if the mean curvature of the initial data set instead sat- 
isfies tr^K < 0, the corresponding existence result applies for the flow with 
speed equal to the reciprocal of H — P, with analogous interpretations of 
the solution in relation to marginally inner trapped surface (MITS) in the 
initial data set. 

The results of this paper are laid out as follows. We begin in Section 2 
with a brief remark on the classical evolution by inverse null mean curva- 
ture, and derive an interior estimate for the null mean curvature of smooth 
solutions. In Section 3 we introduce the level-set formulation of the flow, 
and prove existence of solutions, u e , to the elliptic regularisation problem. 
The translating graphs E| given by ([3]) are then used in Section 4 to study 
the jump regions of the limit, u, of the regularised solutions u £ . In Sec- 
tion 5 we introduce the variational formulation of weak solutions, using the 
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jump region analysis of Section 4 to motivate the choice of definition of weak 
solutions. In Section 6 we introduce the concept of outward optimisation 
to give a geometric characterisation of jump regions of weak solutions, and 
show that the interior of the jump region is foliated by smooth MOTS. The 
main existence result, Theorem [TJ then follows in Section 7, and we discuss 
applications of the flow, including Proposition [2] in Section 8. 

Acknowledgements. I would like to sincerely thank my advisor, Gerhard 
Huisken, for introducing me to this topic and for his insightful and enthusi- 
astic guidance throughout the development of this work. 

1. The smooth flow 

Since the aim of this work is to develop the weak theory for the evolution 
by inverse null mean curvature, we will not provide a classical PDE analysis 
of (*), except to remark that the leading order term of the linearised equation 
is jjj^pyj/S. on the right hand side, where A denotes the Laplace-Beltrami 
operator with respect to the metric g at time t. This is an elliptic operator 
as long as (H + P)~ 2 remains non-singular, so (*) is parabolic so long as the 
null mean curvature of the evolving surface remains strictly positive. 

In Section 3 we construct an explicit, non-compact solution Ef of (*), 
for which we require an upper null mean curvature bound. The objective 
of this section is therefore to derive the interior H + P estimate (J7J) for 
smooth solutions of (*) (which also holds for non-compact solutions). We 
begin by stating the evolution equations for some fundamental quantities. 
Let V be the connection on the initial data set (M, g, K) and let the induced 
connection and second fundamental form on Ef be denoted by D and A = 
{hij} respectively. 

Lemma 3. Smooth solutions of (*) with H + P > satisfy the following 
evolution equations. 
,d 1 \D(H + P)\ 2 1 „ „, - , „ 

d 



in) j f P = {Vutr M K - (V u K)(u, v)) - * D^H + P)K tu . 

iv) — |£t| + / P dV = |Et|, whenever T,q is closed, where V(E) 

dt V(£ t )\V(£ ) 
denotes the volume enclosed by E. 

Proof. The relevant evolution equations satisfied by general flows are recorded 
in [T4"] , except for the evolution of P which satisfies 

^-P= ±tr M K - vV^-Kij - WKiAv* 
dt dt dt dt 

= WTP ^^ mK - & v K){v, u)) - {H + p)2 D t (H + P)K iv . 
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□ 



Combining i) and ii) of Lemma [3] above, we obtain 

d A(H + P) 2\D(H + P)\ 2 \A\ 2 + Ric(u,u) 



() dV ' (H + P) 2 (H + P) 3 H + P 

VjMmK - (V U K) (u, v) 2 A (H + P)K, 



and for the speed function tp : 

8 ^ .,.2/ A-/. , !\ ytl2 



# + P (H + P) 2 

1 



H + P' 



(6) = V 2 (A^+(|^r+ffic(i/, ^)+V i/ tr M iiC-(V^)(^ v))i>+2Diil>K iv ). 
at 

Like in [ID], the supremum o~(x) of radii r for which the interior curvature 
estimate ([7]) below holds is defined as follows. 

Definition 4. Let d x denote the distance to x. Then for any x G M , we 

define u(x) G (0, oo] to be the supremum of radii R such that Br(x) CC M, 
Ric > — looFn^ITR 2 i n Br(x), and there exists a function p G C 2 (Br(x)) 
such that 

p(x) = 0, p > d 2 , on dBft(x), yet \Vp\ < 3d x and V 2 p < 3# on Br(x). 

Lemma 5 (Interior null mean curvature estimate.). Let Sf be a smooth 
solution of (*) on M for < s < t. Then for each x G St and i? < cr(x) 

(7) H(x,t) + P(x,t) <max [ (H + P) R , 



R ( Va 2 + 2nX-a 



where A := 4 (3n + (12 + 3n)||if || c oi2 + n\\K\\ c iR 2 ), a := 12 + 4n||i^|| c oi2 
and (H + P)r is the maximum of H + P on Br, the parabolic boundary of 
E t nB R (x). 

Proof. We wish to construct a subsolution to ([6]). Since 

\A\ 2 > — >- ((H + P) 2 - 2P(H + P)) , 
n n 

and Dip < \Vtp\, P < n||if||(7o and V„P < 7i||2if from ([6]) we obtain 

(8) ^ >^ 2 AV + - - inn , p2 - 2\\K\\c^ 2 

at n WU{n + \)R A 

- n\K\ c ii\) Z - 2|V^| HKllc-oV 2 . 

We allow the evolving surface St to have a smooth boundary <9Sf, and define 
the parabolic boundary of the flow St n Br to be 

B R = B R (x, t) := (B R n S ) x {0} U (U < s < t (B R n <9S S ) x {a}), 
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(H + P) R = (H + P) R (x,t):= sup H(y,s) + P(y,s). 

Consider the function (ft = 4>s{y) := -77 (-R 2 —p(y)), where 

R 

C s := (max(R(H + P) R , A - S, 



for < S <C 1 and p as defined above. Note that A(ft = tr Si (V 2 (ft)-H (V0, ^) . 
Then for y G E t n -Br, we have 

9 ,2A^i (hj ^ 



— — f Al ^ = (V0, ^) - ^tr E( V^ + tf^(V0, v) 
0) = -§(V P> »/> (</> + £ - P<ft 2 ) + 2 §tr SinBfl (V 2 p). 

Since < CxR < 5R < tb, it follows that (ft < ib on Br. In 

~ (# + P) R 

order to obtain a contradiction, let < s < t denote the first time when 
(ip — (f))(y, s) = for y G E s n Br(x). At this point 

^ - ^ 2 a) - 0) < o. 

On the other hand, since (ft < R, it follows from ([8]), Q and the conditions 
on p defined above that at the point (y, s) 

>J i- - 2||K|| c o«^ - nll^Hc^ 2 - 2|V<A| ||K|bo0 + ^(Vp, i/) (2 - P(ft) 

- ^^tr StnBj j(V 2 p) 
-^l2n ~ 2C &(^> + n W K \\coR) ~ C|(3n + 12||K|| c o J R + ?i||ir|| c i J R 2 

+ 3n||ir|| c oi?)J = 0. 

Thus ip > (ft on all of n Br{x). In particular ^(x, t) > (ft(x, t) = C$R, and 
as 5 was arbitrary it follows that ift(x,t) > CqR. □ 

In Section [2] we see that the null mean curvature upper bound given 
by Lemma [5] is the key to existence and regularity, and that this estimate 
continues to hold for weak solutions. On the other hand, the reaction term 
— j/^jp in the evolution ([5]) of the null mean curvature in general leads to 
singularity formation in finite time, analogous to inverse mean curvature 
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flow. We therefore turn to the question of a weak formulation of solutions 
to the evolution by inverse null mean curvature. 



2. Level-set description and elliptic regularisation 

In this section we outline a level-set description of the evolution by inverse 
null mean curvature. This level set formulation allows jumps in a natural 
way, because if u is constant on an open set CI, the level sets "jump" across 
Cl. We use the method of elliptic regularisation as a tool to approximate 
solutions of the degenerate elliptic level set problem by smooth solutions of 
a strictly elliptic equation. Studying the properties of the regularised solu- 
tions helps to guide us towards the optimal formulation for weak solutions 
of (**), which we then define in Section |H 

Level-Set Formulation. The following ansatz lies at the foundation of 
the level-set formulation. We assume that the evolving surfaces are given 
by the level-sets of a scalar function u : M — > K via 

(10) E t := {x : u(x) < t}, S t := dE t . 

Employing the terminology coined by Brian White in [T7], we call u the 
time- of- arrival function for the evolution by null mean curvature. Then 
wherever u is smooth and Vu 7^ 0, the normal vector to St is given by 
Vu 

v = 1 and the degenerate elliptic boundary value problem 

I Vu\ 

u =0, 

describes the evolution of the level-sets of u by inverse null mean curva- 
ture. In this smooth setting, the left hand side is the null mean curvature 
of Y*t and the right hand side is the inverse speed of the family of level-sets. 
Since | Vit| = H + P, the local uniform estimate ([7|) for the null mean cur- 
vature suggests that it is reasonable to expect locally Lipschitz solutions of 
(**). However, in order to interpret (**) as the level-set formulation of the 
classical, expanding flow (*), it is necessary for the zero function be a sub- 
solution barrier for the Dirichlet problem (**). In particular, this suggests 
that it only makes sense to study (**) on initial data sets (M,g,K) satisfy- 
ing ItmK > on M\Eq. We see below that this mean curvature restriction 
is also necessary for the elliptic regularisation problem. 

Elliptic regularisation. In order to solve the degenerate elliptic prob- 
lem (**), we study solutions of the following strictly elliptic equation on the 



EVOLVING HYPERSURFACES BY THEIR INVERSE NULL MEAN CURVATURE 9 



precompact domain Ql := Fl\Eq, 

(*)* 




on dEo, 
on dF L , 

where Fl '■= {v < L} for an appropriate comparison function v defined 
below. In this section we prove existence of a smooth solution of (*) e . 

Rescaling (*) E via u £ := — gives 



and we see that the left hand side corresponds to the null mean curvature 
H £ +P £ of the hypersurface graph(u e ) in the product manifold (M n+1 xl, g), 
where g := g © dz 2 and the given data K is extended to be constant in the 
z-direction. This rescaled equation (*) £ has the geometric interpretation 
that the downward translating graph 

(11) E| := graph ^u £ - t - 

solves (*) smoothly in fix, x 1R. This is equivalent to the statement that the 
function 

(12) U £ (x, z) := u £ (x) — ez, (x,z)gQ l xR, 

solves (**) in x since Ef = {U £ = t}, that is, U £ is the time-of-arrival 
function for the solution Ef. Therefore elliptic regularisation allows one to 
approximate solutions of (**) by smooth, noncompact solutions of (*) one 
dimension higher. 

In fact, (*)g has the further interpretation as Jang's equation ([1]) with a 
gradient regularisation term. In [T2], Jang used equation ([Tj) to generalise 
Geroch's [7] approach to proving the positive mass theorem from the time 
symmetric case to the general case. He noted, however, that the equation 
cannot be solved in general, leaving the question of existence and regular- 
ity of solutions open. The analytical difficulty is the lack of an a priori 
estimate for sup \ w\ due to the presence of the zero order term tr^ (K). In 
[15] . Schoen and Yau bypass this issue using a positive capillarity regularisa- 
tion term which provides a direct sup estimate via the maximum principle. 
In the case of the Dirichlet problem for Jang's equation, appropriate trap- 
ping assumptions must be placed on the boundaries in order to obtain the 
required boundary gradient estimates (see [2] and [6] , or pQ for an overview) . 



In the case of the Dirichlet problem (*) e , we see below that the zero order 
term trM{K) obstructs the existence of a subsolution barrier at the inner 
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boundary. In order to obtain the required boundary gradient estimate at 
this inner boundary, we must impose the ambient mean curvature restric- 
tion mentioned previously, that tr m{K) = g %3 Kij is nonnegative on M\Eq. 
Similarly, it was observed by J. Metzger [13] that restricting to trj^if > in 
the capillarity regularised problem prevents the solution from blowing-up to 
negative infinity over marginally inner trapped surfaces in the initial data 
set. 

A priori estimates and existence for (*) e . As stated above, we 
will use a comparison function v to prescribe the outer boundary dFi of the 
annulus domain £Il for the Dirichlet problem (*) £ . Since M is asymptotically 
flat, outside some compact set O C M we can choose a radial coordinate 
chart such that for an appropriately chosen a > 0, the function v = alogr 
is a smooth subsolution of the following approximating level-set equation 

W^\) +s \ 9 ~^V^ 

for s G [0, 1] in this asymptotic region M\Q. Let 



E- u £ , s := div | ; ,„ ; = ]+8[g 13 - - 



y/\Vu E , s 



2 



To prove existence of solutions to the Dirichlet problem (*) e , we then con- 
sider solutions of the family of approximating equations 

'E E ' s u £:S = mU L , 
(*)e,s < u £tS = on dE , 

u £)S = s{L-2) on dF L , 

for s G [0, 1], where the subsolution v = a log r prescribes the outer boundary 
8Fl = d{v < L} for both the Dirichlet problems (*) e , s and (*) e . We use 
barrier functions at the inner and outer boundaries to derive the following 
interior and boundary gradient estimates. Aside from the supersolution 
barrier at the outer boundary, the following Lemma follows essentially as in 
HQ1 Lemma 3.4]. 

Lemma 6. For every L > 0, there exists e(L) > such that for < e < e{L) 
and s G [0,1], a smooth solution of (*) e s on &l satisfies the following a 
priori estimates: 

(14) u e , s > -e in Q L , u E)S > v + (s - 1)(L - 2) - 2 in F L \F , 

(15) u £ <C(L,\\K\\ Co ) inn L , 

(16) \Vu £tS \ < H + e + n\p\ on 3E , |Vu £iS | < C(L, \\K \\ c o) on dFi, 

(17) \Vu £jS (x)\ < max | Vu £ , x \ + e + C, xgCIl, 

( 18 ) \ u e,s\ c ^(n L ) < C(e,L,n, \\K\\ c o, \\K\\ c i). 
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Proof. Let |A| denote the size of the largest eigenvalue of Kij on £1^. 
la. We construct a subsolution that bridges from Eq to where v starts in 
the asymptotic region, which allows for unrestricted jumps in the compact 
part of the manifold. 

Let Cut(E'o) be the cut locus of Eq in M. We construct a subsolution to 
(*) e on M\(E U Cut (.Eo))- Define G = E , G d := {x : dist(x,E Q ) < d} 
and choose di large enough that G dL 2 Fl- In general, for a surface moving 
in normal direction with speed /, the evolution of the mean curvature is 
given by 

(19) — = -Af-\A\ 2 f-Ric(v,v)f. 

We can therefore estimate the mean curvature of the surfaces dGd via 
f) ff 

— = -\A\ 2 - Ric{v, v) < Ci(L) on dG d \Cut(E ), < d < d L , 
yielding 

H dGd < maxff + + dd < C 2 (L) on dG d \Cut(E ), < d < d L , 

8E 

where H + = max(0, H). Consider the prospective subsolution 

Vl (x) := /(d) = /(dist(x,G )), x G G dL \E , f < 0. 



Since Vt>i = f'v, we have V^-vi = f'(V e M, ej) = f'hij and thus 



2 
ij 

(g ij - i/V)v?,«i = f'H dGd > fc 2 . 



72 

y 

Hence 

y/f^+^^'v! > -\f'\C 2 + + syfFT&gVKa - \f'\\K vv \ - /' 2 - e\ 

In order to obtain an appropriate subsolution, it is in fact necessary to 
restrict to initial data sets (M, g, K) with g 13 > 0, so that we can discard 
the bad term s\/ f' 2 + £ 2 g lJ K^. Then we can use the following barrier 

f(d) := l(-l + e - Ad ) onO <d<d L . 

If we restrict e such that e < e~ AdL , then \ f'\ = ee~ Ad > e 2 and e 2 < \f'\ < 
e. Then taking A := 2(C 2 + |A| + 2) we obtain 

(f 2 + e 2 )(\f'\C 2 + l/'H^I + f' 2 + e 2 ) < 2e 2 (C 2 + \K W \ + 2)|/'| < e 2 f". 

This shows that the function v ljS (x) := — — ^— ( e -(2C* 2 +|A|+2)<i _ ^ 

2(62 + I A| +2) 

is a smooth subsolution for E £,s on G dL \(EoUCut(Eo)^ for sufficiently small 
e. Furthermore, v\ )S is a viscosity subsolution of E £,s on all of G dL \Eo. Since 
ti > v 1 on the boundary, it follows by the maximum principle for viscosity 
solutions [U Thm 3.3] that 

— du dv\ 

(20) u > v\ > — e in fii, and — > — — > — e on 8E . 

ov ov 
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b) We construct a subsolution on Fl\Fq. Assume L > 1 and consider the 
function v 2 := -^-^v-l + (s-l)(L-2). Then E°' s v 2 = E°' s v + y|Vu| > 

on Fl\Fq. Since the domain is compact, for all sufficiently small e we obtain 
E e,s V2 > 0. From (j20|) we have that u > —e in CIl, thus it > v 2 on c^o and 
u = s(L — 2) = i)2 on <9-Fl- It then follows from the maximum principle that 
u > v 2 > v + (s - 1)(L - 2) - 2 in F L \F , thus 

(21) ^ > -G(L) on 8F L . 

ov 

A rescaled version of v 2 provides the required barrier when L < 1. 

c) The zero order term trM(K) prevents constant functions from being su- 
persolutions to (*) £)S , like in inverse mean curvature flow. We therefore con- 
struct a linear supersolution to (*) £ on Fl\ (EoL)Cut(dFi)) , where Cut(9Ff,) 
is the cut locus of 8Fl in Fl- Consider v%{x) := f{d) = /(dist(x, Go)) where 
Go := 9Fl, Gd := {x : dist{x,dFi) < d} and choose do large enough that 

f) FT 

G do D dE Q . From $HJ) we find — > -Gi(L) on dG d \Cut{dF L ), for 
< d < do, yielding 

H 9Gd > - max#_ - Cid > -C 2 {L) on dG d \Cut(dF L ), < d < d , 

8F L 

where H- = — min(H, 0). Setting v^(x) = f(d) := s(L — 2) + ( m+ — ) d, 



where m > is to be chosen, we obtain 



do 



Vf' 2 + eE?>*f(d) < f{C 2 + 2\g^Kij\ + \K VV \ - f) 



Setting m := C 2 + 2\g^ + |A| ensures y 7 /' 2 + eE £ < s f(d) < for all 
sufficiently small e (so that £</'). Then ^3(2;) is a smooth supersolution 
on Gd \(EoUCut(dFL)) . Furthermore, U3 is a viscosity subsolution on all of 
G do \Eo. Since u = f on di^ and u < f on <9i?0) it follows by the maximum 
principle for viscosity solutions that 

(22) u<f<L + md in Q L , 

du 2 

(23) —<C 2 + 2\gVK ij \ + \p\ + — = C(L,\\K\\ c0 ) on 8F L . 
ov do 

d) Choose a smooth function V4 which vanishes on 8Eq such that 

(24) H + + n\\K\\ c o < <H + + e + n\\K\\ c o along dE . 

ov 

Let v be the normal vector to dEo^ and r be the tangent to dEo, which 
satisfies 

Vv 4 



v = A ._ 1 + \/ X 2 — It, for some A > 1. 
|V«4| 
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thus along 8Eq we obtain 

E°' s v 4 < -H dEo + n\\K\\ c o - (H+ + n\\K\\ c o) < 0. 
A 

This implies that E°v A < in the neighbourhood U := {0 < v A < 5}, for 

sufficiently small 5 > 0. Now define the scaled-up function v§ := — , 

1 — v A /o 

for x € U. So ^5 — > oo for v A — > 5, ie on 8U\dEo, and E 0,s v§ = E 0,S V4 : + 
|Vt, 4 | - |V " 4 ' < E^v A < 0. 

(1 - v A /5y 

For e sufficiently small (depending on L and m) we obtain that E £ ' s vs < 
on the set V := {0 < ^5 < L + m do}. From (|22p we have u < L + m do on 
fj^, thus u < U5 on 5V. 

Then by the maximum principle, u < V5 on V, and therefore 

(25) £<^ = ^<2r ++e + n |A| on^o 

ov ov ov 

for sufficiently small e. 



2) The desired interior gradient estimate (|17|) can be obtained from the 
interior estimate for H + P in Lemma Since we can not apply the result 
directly to (*) E ,s (except when s = 1), we instead rework the proof of Lemma 
[5] for the evolution equation 

dF 1 
to obtain the corresponding estimate 

A 



(26) H(x, t) + sP(x, t) < max (H + sP)r, 



y/a 2 + 2n\ 



Here A and a are as defined above, and (H +sP)r is the maximum of H+sP 
on Pji, the parabolic boundary of n Br(x). 

Analogous to (jlip . the downward translating graph := graph ^— — ~\ 

is a smooth solution of (*) s , described by the level-set function U EtS (x, z) := 
u £ ^ s {x) — ez, since tf t ' s = {U £ ^ s = t}. We then relate estimate (i26|) to |Vu £)S | 
via (*)e,s, which asserts that 



{H + sP) tt , s = ^\Vu £iS \ 2 + e 2 
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Now let B := B™ +1 (x,z), be an (n+l)-dimensional ball centered at the point 

(x, z) C M x R. Since S^' s is a translating solution to (* s ), its parabolic 
boundary is just a translation of 8VLl in time. Furthermore, as |Vu e>s | is 
independent of z, applying (f26|) to £^' s n B yields 



Vu e s | 2 + e 2 < sup max 4/ |Vu e s | 2 + £ 2 + C < max |Vu 6S |+e + C7, 
where C := , , ,, A = — 7 is defined in Lemma [5l For e small enough, we 

( Vo +2nA— ol 

obtain from the boundary gradient estimates 

(27) \Vu EiS {x)\ < max H + + 2 + C(L,\\K\\ c o) + C, 

dE DB R (x) 

which leads to the Lipschitz estimate 

\ u e,s\co,i(fi L ) < C(L,n, \\K\\ c i). 

Then by reworking the proof of the Nash-Moser-De Giorgi estimate ([8], 
Thm 13.2), we obtain 

\ u e,s\c^(n L ) < C(e,L,n, \\K\\ c i), 

for some a = a(£li)- This implies a bound on the Holder modulus of 
continuity for the coefficients of E £,s u, so Schauder theory improves this 
estimate to C 2,a 

(28) \u £ ,s\c^(n L ) < C{e,L,n, \\K\\ c i). 

□ 

Lemma 7 (Existence for the regularised problem). A smooth solution of 
(*) £ exists. 

Proof. We first prove there is a solution of (* £ ,s) f° r s = and small e. Let 
u = and rewrite (*) e ,o as 



Me < 
The map 



F(u) := — = div ( — . ^ = | = e in CIl, 

VlW+1 \^/\Vu\ 2 + l) 

u = on dflL- 



F : C^ a (n L ) -> C Q (ft L ) 



is C 1 , and possesses the solution -F(O) = for e = 0. The linearisation of -F 
at u = is 

PF|o = A 9 : C 2 ' a (n L ) -> C a (0 L ). 

The Laplacian on M is an isomorphism, so by the Implicit Function Theorem 
there exists £q > such that has a unique solution for < e < £q. 



We now fix s £ (0,£q) and vary s. Let I be the set of s such that (*) £iS has 
a solution n e>s £ C 2 ' a (f2i). We have shown that / contains 0. We first show 
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that I is open. Let tt be the boundary value map u i— > u\qq. Consider the 
map 

G : C 2 ' a (n L ) x E -> C a (0 L ) x C 2 ' a (dft L ), 

defined by G(io, s) := = ^ £ ' s (w), vr(u;) - s(L - 2) X dF^j , so that 

(*) e , s is equivalent to G s {w) = (0,0). G s (w) is C , and possesses the so- 
lution G°(u ) = (0,0), where uo is the C ' a (Q,L) solution from above. The 
linearisation of G° at uq is the operator T>G\ Q given by 

(29) VG° U0 = ^ A ^ ViV ^+ j : c 2 < a ->• C7°(^ L ) x C 2 ' Q (<9ft L ), 
where 



•yi + ivuoi 2 V i + |v«oiv' v1vwF+T' 

Since P£'2 (^) is a linear elliptic equation with Holder continuous coeffi- 
cients, we can apply Schauder theory (eg [8], Thm 6.14) to deduce that 
2?G* is a bijective map with continuous inverse. It follows from the Im- 
plicit Function Theorem that G maps a neighbourhood of (uo , 0) onto a 
neighbourhood of (0,0). Thus I is relatively open, which completes the 
proof of existence of u e 6 C 2,a (ti) solving (*) e . Smoothness then follows 
from standard Schauder estimates. □ 

In view of the local uniform Lipschitz estimates for u e , by the Arzela 
Ascoli theorem there exist sequences e% — > 0, Lj — > oo, a subsequence Ui and 
a locally Lipschitz function u : M\Eq — > R such that 

(30) Uj — >■ u 

locally uniformly on M\Eq, and by (|27p . « satisfies 

(31) |Vu(x)| < sup H + + C(L,n, \\K\\ c i,R). 

dE nB R (x) 

In the next section we will study the limit of the translating graphs = 
{U e = t}, where the time-of- arrival function U e was defined by (|12|) . By 
setting 

(32) U(x,z) := u(x), 

we obtain that U —> U locally uniformly on (M\Eq) x R, therefore U is the 
time of arrival function of the limit of the smooth flow 1 1— > S|. 



3. The limit of the translating e-graphs Sf. 

Our choice of variational formulation to define weak solutions to (**), 
detailed in the next section, is motivated by: 

1. The variational properties of smooth solutions of (**), 
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2. The limiting behaviour of the family S| of translating solutions of (**) 
in M x R. 

In particular, we show that the sets E t = {u < t}, associated to a smooth 
solution u of (**), minimise the following parametric energy functional 

(33) j£ v {F) := \d*F n A\ - [ \Vu\ - (^ - !/V) K ih 

J FnA 

for each t > 0. That is 

(34) j£„(Et) < j£ v (F), 

for each set F of locally finite perimeter that differs from the set E t on 
a compact subset A of the domain. Here P = trs t if = (g lJ — v l v 3 )Kij, 
where v represents the unit normal Vu/|Vw| to the surfaces T, t = dE t . The 
functional (j33|) . together with the minimisation principle (I34h . generalises 
the variational formulation employed by Huisken and Ilmanen in [10], and 
accordingly allows the evolving surfaces to jump instantaneously over a pos- 
itive volume at plateaus of the time-of-arrival function, u. However, in the 
weak setting, Vit/|Vu| is undefined on plateaus of the locally Lipschitz func- 
tion u, so in order to incorporate the extra P term for this new flow, we 
must define an appropriate notion of normal vector in these jump regions. 
In this section we show that such a vector field can be obtained by taking an 
appropriate limit of the translating graphs Sf. Since the null mean curva- 
ture of these surfaces is uniformly bounded, results of measure theory allow 
us to control them in C 1,a , which leads to a foliation of the interior of the 
jump region {U = to} = {u = to x R}, at jump times to, by hypersurfaces 
satisfying the following result. 

Proposition 8. Let U(x,z) = u{x), where u £ C,'(M\Eq) is the limit 
of the solution u e of (**) e , as in (Sty). Then the interior, ICt , of the jump 
region {u = to} xR = {U = to}, at jump times to is foliated by hypersurfaces 
with local uniform C 1,a estimates, where each such hypersurface is either a 
vertical cylinder or a graph over an open subset of {u = to}. Furthermore, 
each hypersurface bounds a Caccioppoli set that minimises Ju,u in Kt Q , where 
v denotes the C?'° normal vector field to the hypersurface foliation. 

The normal vector held v to this foliation extends — r = 7— , across 

|W| |Vu| 

the jump region ICt in M x R, and this extended vector field helps motivate 
the definition of weak solutions to (**) in the next section. In this context, 
hypersurfaces and sets in M x R will be denoted by the ~ superscript for 
the remainder of this work, unless otherwise stated, and V denotes the 
connection on (M x R, g). 

To prove Proposition we utilise the following compactness result for 
sequences of minimisers of (j33|) . 

Compactness Property 9. Let £1 C M x R , and let Ei C (l be a sequence 
of sets with C]^ boundary such that dE — > dE, locally in C 1,a , with outward 
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unit normal Vi £ (7?'" (TO) to dEi satisfying vi — > v locally uniformly. Let 
Ui S C?,'c(^) satisfy Ui —> U locally uniformly, and assume that for each 
A CC Cl, sup A |VC/j| < C(A) for large i. 

If the sequence Ei minimises Jui,vi on (l, then E minimises Ju,v in 
Proof. We use the inequality 

(35) J UitV . (Eh.) + Jvm (E 2 ) > Jt/^ {Ei U £ 2 ) + (Ei n £ 2 ), 

for an appropriate choice of Caccioppoli sets E\ and E-i such that E\C±Ei 
is precompact. 

We first prove that E minimises Jxj,v on the outside in CI. To this end, 
consider F D E with F\E CC and a suitable compact set GcSJ contain- 
ing F\E. Since the boundary of G is not necessarily Lipschitz continuous, 
we consider a compact set G C Cl with smooth boundary and G C int(G) 
such that 

\d*(FUEi)ndG\ = \d*(Fr)Ei)nd*G\ = \dEiC\dG\ = 0, 

for all i, with traces satisfying f 9 ^ l^^g — V 9 ^ |<i% n+1 — >■ 0. This is possible 

because F U Ei E and Ei E in L 1 1 oc (0\G). Then setting Ft := Ei U 
(F n G) we see that 

\d*F l nn\ = \d*E i n(n\G)\ + \d*(PuE l )nG\+ [ \<p-- P -<p+\. 

JdG U hi 1 

Now, since F% is an appropriate comparison function for Ei, we have 
jS^E^^jg ^Fi), implying 

^(^) < u + 1 9 _ |^ - 4J. 

Now inserting £i = Ei and Ei = F into ()35[) we obtain 

(36) > jg^EiHF) - [ _\<p~ -<p+\. 

JdG 1 ' 

Next we pass to limits. Since the trace term converges to zero, using lower 
semicontinuity we obtain 

JuAF) > J$ V (E). 

The fact that E minimises Ju %v on the inside in Qt amongst compet- 
ing sets F C E satisfying E\F CC can similarly be proven by again 
constructing G and considering the comparison function Fi := Ei n F for 
i » 1 large enough. □ 

To prove Proposition [8] we will also draw upon regularity theory for obsta- 
cle problems of the type (|38p below. In particular, if the set E% := {u < t] 
minimises J u v , then it is almost minimal in the sense that 

(37) \d*E t C\B R \ < \d*FnB R \ +C(n,\\Du\\ oo ,\\K\\ c0 )ir +1 , 
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for EtAF CC Br. This means we can apply partial regularity results of 
geometric measure theory to obtain higher regularity for the level-sets T,t = 
dEt. Specifically, we consider the following C 1,a result (see for example 
[16]), as quoted in [TO] , 

Regularity Theorem 10. Let f be a bounded measurable function on a 
domain f2 with smooth metric g and dimension n + 1 < 8. Suppose E 
contains an open set A and minimises the functional 

(38) \d*F\ + [ f 

Jf 

with respect to competitors F such that FDA, and FAE CC S7. If dA is 
C 1 ' 01 , < a < 1/2, then dE is a C 1,a submanifold of Q with C 1,a estimates 
depending only on the distance to dfl, ess sup\f\, C 1,a bounds for dA, and 
C 1 bounds (including positive lower bounds) for the metric g. When n + 1 > 
8, this remains true away from a closed singular set Z of dimension at most 
n — 7 that is disjoint from A. 



Proof of Proposition We break up the proof into the following Lemmata. 
Lemma 11. The level-sets Ef = {U e = t} are locally uniformly bounded in 

Proof. Since Ef = {U e = t} is a smooth solution of (*) on (M\Eq) x R, with 
smooth normal vector field u £ = ^f-\ i the functional Jjj e ,u e is well defined 

for sets F C M x M of locally finite perimeter. Using u e as a calibration 
and applying the divergence theorem exactly as in the proof of Smooth Flow 
Lemma [T71 shows that E\ := {U £ < t} minimises Ju e ,v e hi ^ := E^\E e a for 
a < t < b. 

Now consider x = (x,x') e (M\E ) x R, and d = dist(x,9£'o x R) = 
dist(x, 8Eq). Take L' large enough that B^jix) C Fy . Then for e < e' = 
e(L'), (|27p provides a uniform bound for |Vn e | (and thus also |Vf7 e |+Pg. e ) on 

Bj IxR (x). It then follows from Theorem [TO] that the surfaces Ef n -Bj +1 (x) 
are uniformly bounded in C ,a in e and t. □ 

Lemma 12. Let fCt denote the interior of the jump region {U = to}, at 
a jump time to- Then each point Xq = (xq,zq) G /Q lies in a surface 
Ex C K-t that is locally uniformly bounded in C 1,a , and is either a vertical 
cylinder or a graph over an open subset of /Q . 

Proof. The sought-after surfaces are constructed using a pointwise approach 
similar to that used by Heidusch [9] to prove local uniform C 1, regularity 
estimates for the level-sets of the weak solution to inverse mean curvature 
flow. In particular, we fix a target point 

(39) Xq = (xo, zq) £ /Ct (J , 

and construct a surface containing that point. 
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Given the convergent sequence £j — > that produces the limit u of the 
elliptic regularised solution u e as in ()30|) . we consider the corresponding se- 
quence of times, ij, at which the surfaces = graph ( ^ — ) pass through 

the target point Xq. This is possible because the translating graphs Ti\ for 
— oo < i < oo foliate x R, thus for every i there is a unique U such that 

In order to write each surface SJ. locally as a graph over its tangent space 
Tx XJ., we use the exponential map to work locally in normal coordinate 
charts on small Euclidean balls B n+2 . In particular, let l{Xq) be the injec- 
tivity radius of Xo in M\Eq x R, and set 

(40) d = d(X ) = min(i(Xo),dist(Xo,^ to )). 

By Corollary 1111 there exists £o > such that for all t and £ < £q, the surface 
pieces n-B^ xR (Xo) are uniformly C 1,a bounded in t and e. Now consider 
the exponential map 

(41) exp Xo = (exp X0 ,id R ) : T Xo (M xl)n £^ +2 (0,z ) -> xR (X ), 
and set 

(42) % = exp-^S^ n Bf xR (X )) C T X() (M x R). 

In the R-direction the exponential map is just the identity, thus each surface 
Yj\. translates downwards in exactly the same manner as SJ.. Furthermore, 
the surfaces Sf. are uniformly C 1,a bounded in t and e. 

Then there exists R > 0, depending only on the locally uniform C 1,a 
bound, such that 5™ +2 (X ) C t\. and thus the surface pieces S|.n^ +2 (Xo) 
possess uniform C 1 '" bounds. Here Xo = (xo, zq) = exp~ 1 (Xo) is our target 
point. 

The corresponding normals v^Xq) to T\. D -B^, +2 (Xo) create a sequence, 
a subsequence of which converges uniformly to a vector v{Xq). The normal 
space to v(Xq) defines a hyperplane T containing Xo. Then by taking i 3> 1 
large enough, we can write the converging surfaces ni?^, +2 (Xo) as graphs 
of C 1 '" functions Wi over T. By reducing R, and taking i 3> 1 large enough, 
we can then write each Sj. locally as the graph of Wi over T D B 7 ^~ 1 {xq). 

By Arzela-Ascoli, there exists a further subsequence Wi j and a C 1 function 
w:fn B™ +1 (x ) -> R such that 

(43) ^. -> u) in C^f n ££ +1 (£ )). 

Here w is locally the graph of a surface around Xo, and T = T%llx - 
Since the C 1,a bounds on Wi were independent of i, it follows that w G 
C 1 ' ^ n £™ +1 (f )), with the same uniform C 1,a bounds as Wi. Thus 
exp g (E^- o ) := Ex n-B^ xR (Xo) is uniformly C 1 '" bounded. By successively 
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taking subsequences, the T\, converge to a complete hyper surf ace that we 
will henceforth denote by £x > since it coincides with Sjf near Xq. 

Now Xq G {U = to} where, by hypothesis, to := lim^oo U is a jump 
time. In order to argue that T,x is contained in the set {U = to}, we note 
that it is a consequence of the above construction that any y G Sx is the 
limit of a sequence yi G S|.. The local uniform convergence Ui — > u implies 
that iii —7- u uniformly on £>^ +1 (0). Thus the uniform convergence U — >■ U 
on B^ x ^(Xq), together with the fact that lim^oo yt = y then implies that 
U (y) = to, since 

Pi{Vi) ~ U(y)\ < \Ui(yi) - Ui(y)\ + |t%) - U(y)\ -> 0, 

and lim U l (yi) = lim i, = to, thus Sx C {C7 = to}- 



This approach enables one to choose any point Xq in the jump region 
JCt and construct the corresponding surface Sx containing Xo. Since each 
Sx is the limit of the graphs S|. with local uniform C 1,Q bounds, it is clear 
that each Xq is either a vertical cylinder or a graph over an open subset of 
JQq n M. Therefore, let flc denote the open region in JCt n M where |Vuj| 
converges locally uniformly to a finite limit, and let flc denote the region 
where \Vv,i\ converges to infinity. Then the translating nature of Sf together 
with the above construction dictates that the EJ. converge to a graph Sx 
over Qq, lying in a stack 

(44) {t x J = S Xo + «e n+2 , a G R, 

of vertical translates of Sx - To see this, note that 

~ i ■ f u i U-\ 

X = (xo, zo) G = graph — 7 - 3 - -> graph(w) = Sx , 

V-"-. - / 

implies X Q := (xo, zq + a) G ,_ a£jj . , where 
~i- / u i tj . — a£j . \ 

s 4— ae« = graph ~^ gF>Ph(w) + ae z := T, Xa , 

3 \ "', £i 3 / 

where w := exp q (w). Therefore x IR is bounded by vertical cylinders, and 
filled by the stacks produced by the family of vertical translations of 

each graph T>x - □ 

The possibility of two surfaces Ep x and Ep 2 from Lemma 1121 touching 
tangentially at one point P, such that the outward unit normals agree at P 
and £ p 1 lies outside X p 2 (in the direction of the outward normal near P) is 
ruled out by the strong maximum principle. Furthermore, the intersection 
of two surfaces in the limit is ruled out by the translation invariance of the 
surfaces E| and their local uniform C ,a bounds. 

We now argue that we can construct a "normal" vector field v on JCt 
using the surfaces from the proof of Lemma [T2l Since the limit surfaces 
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are vertical cylinders or stacks of translation invariant graphs, the normal 
vector field v to the family of surfaces in /Q is translation invariant, and 
we need only show that we can construct Ex for each Xq C JCt Pi M. 
Therefore, choose a dense set of points in JCt n M. This corresponds to a 
countable set of points {pi}, and for each such pi E /Q n M, we consider 
the convergent subsequence e, such that E|. converges to the hypersurface 
T,p i in JCt Q , where Pi := (pi,0). Then by taking a diagonal subsequence e^, 
we obtain local convergence of Ej* to Ep. for every point pi in the dense set. 

Now consider a point pq E $7^ such that po is n °t in {Pi}- We wish to 
argue that we obtain local convergence to Ep via the convergent sequence 
£i m . There exists a point pi in the dense subset such that dist(pi,po) < d/W. 
Let 

(45) dp := min(t(P),dist(P,dX: to )). 

By Corollary II 11 the surfaces Ef are uniformly bounded in B^ IxK (Pi). Then 

Pi 

since Pm /io(Po) C Pm p . (Pi), the surfaces Ef C\B d /io(Po) possess the same 
uniform C 1,a bounds and we can take a convergent subsequence of £$„ such 
that we obtain convergence to a limit surface Ep in Prf p ./io(Po)- Therefore 
this approach constructs a complete graph through each point xq G Og, and 
we obtain the vector field P in all of 

Then given the uniform C 0,a normal vector field v of the hypersurfaces 
constructed through the dense set of points {pi}, we can extend the vector 
field v to any points that have been missed in Vic- Then translating u in 
the e n +2 direction, we obtain a normal vector field on the entire jump region 
/Ct . For the remainder of this work, let v denote this translation invariant 
normal vector field to the surfaces Ex foliating Kt . 

Lemma 13. Let u denote the normal vector field to the surfaces foliating 
the jump region JCt , as above. Then each surface Ex in the jump region 
bounds a Caccioppoli set that minimises Ju,v i n K-t ■ 

Proof. Consider the Caccioppoli set E that is bounded by the limit hyper- 
surface Ejfo ) such that v is the outward unit normal of the relative boundary 
dE(~]iCt . The sets E\. minimize the functional J\j uVi in tt to , where z/, = . 
Passing these sets to limits as in the proof of Lemma [12] to obtain the limit 
surface E x > Theorem [9] then says that E minimises Ju^ in )Q . □ 

Collecting the above results, we obtain a family of Ch" hypersurfaces 
foliating x an d by extending the family of cylindrical hypersurfaces 
in x R to any missed points in Qo, we obtain a foliation of the entire 
interior region ICt . At each point Xq = (xo,to), the corresponding leaf of 
the foliation passing through Xq is constructed by taking the limit of the 
Ef locally around Xq, as in Lemma [T2l This completes the proof. □ 
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4. Variational formulation of weak solutions. 

By freezing |Vit| —tr^K and treating it as a bulk term, one may interpret 
(**) as the Euler-Lagrange equation of the functional 

(46) J UtV {v) := J \Vv\ + v (|V«| - {g ij - i/V) K i3 ) dx. 

For a smooth family of solutions of (*), we will see below that the corre- 
sponding time-of-arrival function u defined by (llOj) satisfies 

(47) J u>1/ (u) < JuA v )> 

among competing locally Lipschitz functions v, that differ from u on a com- 
pact subset of M\Eq. The relationship between the variational formulation 
(|47p and the functional (|33|) is then given by the following Lemma. 

Lemma 14. Let u be a locally Lipschitz function in the open set O, and v 
a measurable vector field on TfL Then u satisfies j^7| ) on 0, if and only if 
for each t, E t := {u < t} minimizes [33\) in fi. 

Proof. This follows exactly as in |1CH Lemma 1.1], with | Vu| — (g* J — u l v 3 )Kij 
replacing the bulk term |Vu|. 

□ 



This equivalence between the two variational formulations also extends 
to the initial value problem 

(48) 



u GC,' (M), v a measurable vector field on T(M\Eq), 



E = {u < 0}, and u satisfies (JUD in M\E . 

To see this equivalence, let Et be a nested family of open sets in M, closed 
under ascending union, and define u as in the statement of Lemma [TH by the 
characterisation Et = {u < t}. Then using Lemma [TJ] and approximating 
up to the boundary, we see that (08]) is equivalent to 

u G C lo ' c M, v a measurable vector field on T(M\Eq) 
and Et minimises J u v in M\Eq for each t > 0. 



Lastly, by approximating s \ t, we see that (1481) and (|49l) are equivalent to 
it G Ci°' c (M), i/ a measurable vector field on T(M\Eq) 

(50) 

and {u < t} minimises J UyU in M\Eq for each t > 0. 

We now present the precise definition of weak solutions to (**). In the 
previous section we highlighted the need to define the normal vector field v 
in jump regions in order to incorporate the P = (g 13 — v l v 3 )Kij term into 
a variational formulation of weak solutions to (**). We showed that taking 
an appropriate limit of the smooth translating solutions Ef = {U e = t] 
of (*), provides a constructive method of foliating the interior JCt a of the 
jump region {U = to} = {u = to} xM, one dimension higher, by C lo '" 
hypersurfaces Sx inMxl with uniform C^'° unit normal vector field v. 
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Each such hypersurface Sx in the foliation is either (part of) a vertical 
cylinder, or is a smooth graph over an open subset of /Q , in the stack 

(51) t + ae n+2 , Q£l, 

of vertical translates of S. The normal vector field v to each vertical cylinder 
is perpendicular to the z-direction, and could therefore be projected to M 
without loss of information. However, in the case of the graphical hypersur- 
faces (|5ip , information would be lost if one were to define the vector field v 
in (|33p to be the projection of P to TM. 

This motivates the choice to formulate the weak solution to (**) one di- 
mension higher, in terms of a translation invariant function U(x, z) = u(x) G 
C \ol( M x K )> and a translation invariant vector field P G C^(T((M\E ) x 
R)) that extends V£//|VJ7| across the jump region. One then considers the 
analogously defined functionals Ju,v to (|3B| and ([33]) for such pairs (U, P) in 
M x R, and we remark that Lemma [TH and the initial value problem equiva- 
lences (|48p -(j50l) hold in M x R (for general U and P that are not necessarily 
translation invariant, like we will demand for the weak solution of (**)). 

In Lemma [13] we showed that each of the surfaces Xx foliating the jump 
region ZQ bounds a Caccioppoli set that minimises Ju,i> in the jump region 
lCt - Together with Lemma [HI this motivates the restriction in Definition 
[TBI below that at each point X G (M\Eq) x R, P(X) be the normal vector 
to a C 1,a hypersurface that bounds a Caccioppoli set minimising Jjjp in 
(M\E ) x R. 

Definition 15. Let Eq C M be a precompact, open set with C 2 boundary 
So = OEq. We call the pair (U,P) a weak solution of (**) with initial 
condition E if U G C^(M x R) and v G C^(T(M\E Q ) x R) safe/y 
(i) U is translation invariant in the vertical direction. In particular, there 

exists a locally Lipschitz function u : M — > R such that U(x, z) = u(x) 

and u satisfies 

■ u{x) > Vx G M\E , 

■ u\ dE = 0, u{x) < Vx G Eq, 

■ u(x) — > +oo as dist(x,Eo) — > oo. 

(m) T/ie sei ^ = {[/ < t} minimises Ju t p in (M\Eq) xR for each t > 0. v4i 
jump times to, each point Xq = (xq, zq) in the interior ICt of the jump 
region {U = to} lies in the boundary dEx Q G of a Caccioppoli set 

Ex that minimises Jjj,d in K-t ■ 
(Hi) v is a translation invariant, unit vector field such that 

■ v{X + ae z ) = v(X) MX G (M\E ) x R, a G R, 

• P(X) is the normal vector to dEt at each point X G dEt, 

■ P(X) is the normal vector to dEx at each point X G dEx , 
at jump times tQ. 



24 KRISTEN MOORE 

Remarks 1. Unlike in the weak formulation of inverse mean curvature 
flow, which asks only that Et = {u < t] minimise J U)V for each t > 0, we 
require the variational principle (|34p for Jjj t p to be satisfied everywhere, in 
particular in the interior of the jump region. 

2. By Lemma [T4l any weak solution (U(x,z) := u(x),i/) of (**) satisfies 
([38]) on (M\E~q) x R. Furthermore, we find that (u,um '■= v\ TM ) satisfies 
g7D in M\E . 



Lemma 16. Let (U(x,z) := u{x),v) be a weak solution of (**) with initial 
condition Eq. Then the pair (u,,i>m) satisfies {J^lty on M\Eq, and Et = 
{u < t] minimises J u ,v M for each t > 0, where vm '•— v\ TM - 

Proof. Since the tensor K is extended trivially in the z-direction, we find 
(52) (g ij - pV) Ka = (gV - i/ M v\^ K ij} 



where vm '■= u \tm~ ^ u - v m := \^ u \ ~ \9 % ^ ~ ^m^m) denote the bulk 
term of J u ,v M - Let v be a locally Lipschitz function such that {v / u} C 
A CC M\Eq. Let <t>(z) be a cutoff function such that: 

\4>z\ < 2, (j) = 1 on [0, s] and = on K\(-l, s + 1). 

Then V(x,z) := (j)(z)v(x) + (1 — 4>{z))u(x) is an appropriate comparison 
function for U, and letting A := A x [— 1, s + 1], we obtain from (|47p 

[ \Vu\ +uB u>UM = [ | W| + U (|W| - - PV) 

< [ | W| + V (I W| - - PV) Ky) 

< / 4>{\Vv\+vB u>VM ) + {l-4>){\Vu\ + uB UtVM ) 

J A 



+ k& 2 h - u . 



This implies 

s • Juu M ( u ) = s \^ u \ +uB Ufl 
J A 

< J>(\Vu\+uB u>UM ) 



■MM 



V — U 



< f <j)(\Vv\+vB lhUM ) + 
J A 

< (s + 2) / \Vv\+vB u>UM + / 

JAx(\-l 

<(s + 2)J u A UM (v)+A [ 
J A 



([-1,0]C[1,2]) 



v — u 



\v — u\ 
I A 
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Dividing by s and passing s — > oo proves that the pair (u, vm) satisfies ([37]). 
Lemma [TH then implies that the sets Et := {u < t] minimise J u ,v M for each 
t > 0. □ 

We now state some further properties of weak solutions of (**). We begin 
by showing that smooth solutions of the flow (*) are weak solutions in the 
domain they foliate. This follows as in [10\ Lemma 2.3]. 

Smooth Flow Lemma 17. Let (£t) a <t<6 be a smooth solution of (*) on 
M . Let U = t on x R, U < a in the region bounded by T, a x R, and 
Et := {U < t}. Then Et minimises Ju,u in Eb\E a for a < t < b, where v is 
the smooth normal to the vertical cylinder St x K, given by v = (z/£ t , 0) = 
W 

jw|' 

vu 

Proof. We use the smooth normal v = jy^yj" as a calibration and apply the 

divergence theorem to relate Ju,u{Et) to Ju,d{F) for a competing set F of 
finite perimeter with FAE t CC 0. Let B U)V := |W| - (jfi - u^^Kij 
denote the bulk energy term in Ju,u- 

JuAEt) =\dE t \ - [ B UiP dx = [ v Q&f ■ vdU n ~ x - [ B v . p dx 

JE t JdE t JE t 

v d . ■vdU n ' 1 + I v Q& -vdU n - 1 - f B u>D dx 

dEtHF Et JdE t \F ' JEt 



I 



v d ,p-vdH n l + I Bjj,&dx— / Bjj^dx 

d*FnE t J Et\F JE t 



+ / v^p-vdW 1 - 1 - \ B UtP dx 

Jd*F\E t J F\E t 

= [ v d *r ■ vdU n ~ l - [ B v . p dx < \d*F\ - [ B UyV dx = JuA F )- 

Jd*F JF JF 

□ 

Weak Mean Curvature. In view of the local C 1,a estimates given by 
Regularity Theorem [TO] we can consider the weak mean curvature of the 
surfaces S t = d{U < t}. 

Let N be a C 1 hypersurface in M x BL Then a locally integrable function 
H on N is called the weak mean curvature provided 

(53) f div^Xdfi = [ Hv Xd^i, G C C °°(T(M x R)). 

Jn Jn 

Lemma 18. Let Et := {U < i] minimise Jjj^ in A := Efj\E a , for U G 
C^l(A). Then the surfaces S t = dE t have weak mean curvature H satisfying 
H = |Vf7| —P for a.e. x G £t and a.e. t G (a, b), where P = (g 13 —D l i> :, )Kij . 
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Proof. Let X be a compactly supported vector field denned on M, and 
(<3> s )_ £<s<£ the flow of diffeomorphisms generated by X with <3?o = id,M- 
For minimisers of Ju, v , we use the area formula, the dominated convergence 
theorem and the co-area formula in the form 



p roc r 

/ \Vf\dx= dt 

JW l + 2 J-oo J{f=t} 



to obtain 
d 
ds 



d 
ds 
d 
ds 



J U>V {U o $-1) 

[ \V(U o cD;i)| + (u o (|W| - - dx) 
Jw / 

det d$ s (x)\dH n (x)dt 



s=0 



s =o \J —do Js t nvy 



- / W • X (| VET] - - y^Kij) dx 
Jw 

/CO /■ /* 
/ div^ Xd?T^ - / v ■ X\VU\ (\VU\ - (g ij - ^u^Kij) dx, 
-do Jttnw ' •/vk 

W - — 

since ^ = —= — - when VJ7 7^ 0, and ^ VZ7 = when S/U = 0. Then by the 
|W| 

co-area formula, we obtain 

0=/" / (div g X + (P- |W|)z>) -Xdn n+l dt. 
J-00 Jt t nw ^ * ' 

Lebesgue differentiation and comparison with (|53p yields the result. □ 

Exactly as in the proof of [101 Theorem 2.1], we also obtain the following 
compactness theorem for the time-of-arrival function. 

Compactness Property 19. Let Ui G Cf^flj) and i>i G C?J?(TQi) be a 
sequence of solutions of J^7| ) on open sets Cli C M x R, such that 

(54) f)j ->■ fTj — > £/, ->■ z?, 

locally uniformly, and such that for each j4 CC O, sup^ |VJ7j| < C(j4), /or 
iar^e i, where C(A) is independent ofi. Then (U, i>) solves J^7[ ) onVt. In the 
special case where {Ui,Vi) is a sequence of weak solutions of (**) satisfying 
Definition \15\ then the limit (U,v) is a weak solution of (**). 

5. Geometric characterisation of jump regions 

In this section we introduce the concept of outward optimisation in order 
to give a geometric characterisation of the criterion selecting jump times. 
Since weak solutions (U(x, z) = u(x), v) of (**) are translation invariant and 
the level sets of U are vertical cylinders, this characterisation follows from 
the parametric variational formulation (j34|) for {u,vm '■= ^\tm)- 
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Let SI be an open set in M. Then we call a set E outward optimising {in 
SI) with respect to u, if E minimises 'area plus bulk energy P' on the outside 
in SI. That is, if 

(55) \d*E C\A\< \d*F n A\ + [ (g ij - i/i; J }l< ir 

Jf\e 

for any F containing E such that F\E CC SI, and any compact set A con- 
taining F\E. Here v is a measurable vector field on F\E. The set E is 
then called strictly outward optimising (in S7) if equality in (I55p implies that 

Fnn = fino a.e. 

We use this concept to define the strictly outward optimising hull of a 
measurable set E C Q. Specifically, we define E' = E'q to be the intersec- 
tion of the Lebesgue points of all the strictly outward optimising sets in S7 
that contain E. We call E' the strictly outward optimising hull of E (in SI). 
Up to a set of measure zero, E' may be realised by a countable intersection, 
so E' is strictly outward optimising, and open. 

We then obtain the following interpretation of the variational formulation. 

Outward Optimising Property 20. Suppose that (U(x,z) := u(x),v) is 
a weak solution of (**) with initial condition Eq, and that M has no compact 
components. Then: 

(i) For t > 0, Et is outward optimising in M with respect to vm '■= v\ TM - 
(ii) For t > 0, Ef is strictly outward optimising in M with respect to vm- 
(in) For t > 0, E[ = E^ , provided Ef is precompact. 

(iv) For t > 0, \dEt\ = \dEf\ + f E +^(g l: > — v^v^Kij, provided that 

Ef is precompact. This extends to t = precisely if Eq is outward 
optimising. 

Furthermore, for general (U, v) satisfying J^£p in M x K ; the analogous 
statements hold on compact sets Q C M x M. with 

(56) \d*E n A\ < \d*F n A\ + / (g ij - v l v^)K ih 

Jf\e 

replacing \55}) in the definition of outward optimising. 

To prove Outward Optimising Property [20] we will need the following 
Lemma. 

Lemma 21. Let (U,i>) satisfy ft47\ ) on S7. Then U has no strict local maxima 
or minima on £1. 

Proof. First assume that U possesses a strict local maximum so that there 
is a connected, precompact component E of {U > t} for some t. Define the 
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Lipschitz function V k by 
(57) V k := 



k on E k := E k n E, 
U on n\E k , 



for < k < supU and E k := {U > k}. Then (j47j) and Holder's inequality 
yield 

71-1 

(58) [ \WU\(l + U-k)< [ (U-k)C <C ( f (U-k)^) " \E k \$, 

JE k JE k \JE k / 

where Co = (n + 1)|A| and |A| is the size of the largest eigenvalue of K on 
£1. Then using the Sobolev inequality on the left hand side we obtain 

(59) / \VU\(l + U-k)> f |W|= / \V(U-k)\ > ( f (U-k)«- 

JEk JE k JEk \J Ek 

Combining (|58p and (|59p we find 1 < Co|-Efc| 1//n , which leads to a contradic- 
tion since \E k \ can be made arbitrarily small by choosing k close to snpU. 

E 

Now assume that U possesses a strict local minimum and let E be a 
connected, precompact component of {U < t} for some t, and again consider 
the function Vk defined by (I57D . where this time k > inf U and E k := {U < 

E 

k}. Then as above, (|47p and Holder's inequality yield 

n-l 

(60) / |W|(l + ?7-fc)<C7o( / (U-k)^) " \E k fi, 

JE k \JE k J 

and by restricting to k small enough that 1 + U — k > | on E k , we obtain 



n-l 



(61) / \vu\{l + U-k)>Uf_ (U-k)^A 

J E k \J E k J 

Combining ([60|) and (jBT]) we find 1/2 < Co | E k \ , which leads to a contra- 
diction since \E k \ can be made arbitrarily small by choosing k close to inf U. 

E 

In the case where (U(x,z) := u(x),u) is a weak solution of (**), repeating 
the above calculation for u on M, using (|16p . yields the desired result. □ 



Proof of Outward Optimising Property \2(h Refer to [101 Minimising Hull 
Property 1.4]. 

(i) This follows immediately from Lemma [T6l 
{%%) From (|50p we obtain for suitable A 

(62) \d*Ef r\A\< \d*F C)A\+ J (g ij - lAjvi ,)!<,, - \Vu\dx, 

for any t > 0, any F with FAE^~ CC M\E^~ , proving that Ef is outward 
optimising. 
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To prove strictly minimising, suppose F contains Ef and 

\dE+ n A\ - \d*F n A\ = I tfi - u^^Ktj. 

Jf\e+ 

Then by (162H . Vu = a.e. on F\E^~ . Since F is also outward optimising, 
and the Lebesgue points of an outward optimising set form an open set in 
M, by a measure zero modification we may assume F is open. Then u is 
constant on each connected component of the open set F\{u < t}. Since 
M has no compact components, Lemma [21] (i) means that no connected 
component of F can have closure disjoint from E^ , therefore u = t on F\Et 
and F C Ef . This proves that Ef is strictly outward optimising. 
(Hi) It is clear from part (ii) and the definition of E' t that E' t C Ef. Assume 
Ef is precompact. Then if 

\E' t nA\ = \E+ n A\ + / - ^h)K tJ , 

Je+\e[ 

strict outward optimisation implies that E' t = Ef . Otherwise 

\dE' t n A\ < \dE+ n A\ + I tf* - u^ M )K ih 

Je+\e[ 



contradicting (I62j) . 

(iv) In view of (i), we can use Ef as a competitor to obtain 

(63) \8Et HA\< \8E+ n A\ + J (g ij - v^v^K^dx, 

Et\Et 

for t > 0, and for t = if Eq happens to be outward optimising itself. Since 
Ef is precompact, strict inequality in f|63|) would contradict (Hi), implying 
equality in (|6"3"j) . which proves (iv). 

The proof for general (U, v) satisfying (|48p inflcMxl follows exactly as 
above. □ 

Outward Optimising Lemma [201 implies that dEt satisfies the obstacle prob- 
lem minimising "area plus bulk energy P", with Ef as the obstacle. This 
leads to a heuristic interpretation of the minimisation principle (|47p . Namely, 
as long as Et remains strictly outward optimising, it evolves by inverse null 
mean curvature, and when this condition is violated, Et jumps to E' t and 
continues. This implies that the null mean curvature is nonnegative on the 
weak solution after time zero. Furthermore, part (iv) of Lemma [20l implies 
that the monotonicity property 

(64) / P=|£,| 

E t \E 

derived in Lemma [3j also holds in the weak setting, as long as £4 remains 
compact. 
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The outward optimising property also implies a stronger result for the 
surfaces foliating the jump region in Proposition \2E[ namely we see that 
each Sx is a smooth MOTS in K,t . 

Proposition 22. Each surface £x from Proposition in the foliation of 
the interior ICt of the jump region in M xM is a smooth MOTS. 

To prove Proposition [22l we require the following Lemma. 

Lemma 23. 

(65) |Wi|->0 inL} oc {ic t0 ). 

Proof. Recall d defined by (jlOj) . consider a target point Xq = (xq,zq) such 
that zq > 2d + 1 and select a cutoff function 6 G C 2 (R) such that 6 > and 
spt 6 C [zo — 2d, 2o + 2d]. Then let To = zq — 2d — 1, fix an arbitrary time 
T > To, and consider Tq < t <T and L>T + 3 + zo + 2d. 
We wish to show that 

liminf / \D(H + P)\ 2 < oo. 

Jti nB™ xR (x ) 



To this end, we calculate 

\2 



(66) ± [j(z)(H + Pf 

= L t 2 ^ H + P) ¥t {H + p) + {h + p)2 Tz-hTp- + mH + P) 

(H + P)A (jf^p) + \A\ 2 + Ric(u £ , u £ ) - V Ue P 




2Di(H + P) 




+ (H + P)-^ ■ u e + 6H(H + P) 



-2\A\ 2 -2Ric(u £ ,u £ ) + H(H + P) 

DAH + P) \ 6 D(H + P) , ,d6 
+ 2V„,P - ^L—lK^j - 2± . -L-^ + (H + p> J: . „. 

In view of the sup estimates (fb£|) and (f22j) for u £ , there is R(T) > depend- 
ing only on the subsolution v and K j,- such that 

S| n (M x S pt0) C S(T) := (B R(T) \E Q ) x [z - 2d, z + 2d], T < t < T. 

The Outward Optimising Property 1201 applied to E\ compared to the per- 
turbation Ef U S(T), then provides the area estimate 

(67) |Ef n (M x S pt0)| < C(T) + / P < C(T, ||i^|| c o), T <t<T. 

Js(T)\ei 
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Together with the interior estimate ([7]) , and the boundary gradient estimates 
for u £ , this shows 

\H + P\ < C(T, \\K\\ c i) on ff t n (M x spt</>), T <t<T. 

It follows that 

/ ( j ) \H\(H+P)+0(H+P) 2 +\(H+P)V0-iJ £ \ < C(T, \\K\\ci), T <t<T. 
We estimate the D(p and K^ v terms via 



2D4, D(H + P) 



H + P 
H + P 



LD0p |D(iJ + P)| 2 ^ |D(ff + P)| 2 

" 2 (H + P) 2 " 2 (# + P) 2 



2 (H + P) 2 ' 



Thus (1661) becomes 

12 



and integrating gives 



W ft m±ff< C{ T Amcl) , 

JT JT, e t n(Mx[z -2d,z +2d]) K H +^) 

using a <j> such that <j) = 1 on [zq — 2d, zq + 2d] . 

Applying Fatou's Lemma, for any sequence £» — > we obtain 

' \D(H + P ^ 2 

't\C\(Mx[z -2d,z +2d\) (H + P) 



(70) liminf / 1 | ^ < oo, a.e. t > T . 



Now consider the subsequence £j. — > from (|43|) such that E*V — >■ So 

in C^T n S^ +1 (X )), where T = T Xa t Xo . We write t = ij henceforth. 
Since (|70p only holds for a.e. t > To, it will take more work to argue that 
liminf (kj r mxr, y n |-D(P + P)\ 2 < oo. To this end, we pick a sequence t{ 

such that i i — > to + 5 for some \5\ > 0, \U — U\ < e^d and 
(71) liminf / \D(H + P)\ 2 < oo. 

El n(Mx[z -2d,z +2d]) 

Define Zi := ^ °^ — — and 5i := Zi — zq. Then the fact that is just a 

£j £i * ! 
translation of by <5, in the z-direction implies that 

j \D(H + P)\ 2 = J \D(H + P)\ 2 , 

Ej nflf xR (x ) fij nB d MxM (x 0!2o +^) 
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for each i. Furthermore, the condition |tj —U\ < e^d implies that \bi\ = \z 



zq\ < d, which ensures that E| n Bf XK (x , z + Si) CMx [z - 2d, z + 2d], 
and thus from ([71 p we obtain that 

liminf/ \D{H + P)\ 2 <( \D(H + P)\ 2 < oo, 

Jtl nB^'^ixo.zo+Si) Jt i l n(Mx[zo-2d,z +2d]) 

from which our desired estimate follows 

(72) liminf/" |£>(# + P)\ 2 < oo. 

As in the proof of Lemma (|12p . the converging surfaces SJ. can be written 
locally, via the exponential map, as graphs of Cj^" functions Wi over the 
hyperplane T. This local C 1,a convergence of the hyper surf aces, together 
with the first variation of area formula and the Riesz Representation Theo- 
rem then implies that exists weakly as a locally L 1 function, with the 

x o 

weak convergence 

(73) / H ti u & X^f H t v t -X, XeC°(T(M\E x 

H X 

Then by (fT2j) and Rellich's theorem there exists a subsequence (again de- 
noted by i) such that 

in I 2 ^^ nn+h 

^ x o 

Now the level-sets = {U = t{\ smoothly solve (*) in Qi x R, thus 



(74) (H + P) t ^(H + P) t mL 2 (TnB n R ^(X ))). 

h = t>i} smoothly so 
(H + P) & = \VUi\, 



and 

|2 / fxj , D \2 



(75) / \VUi\ z = / (H + py^ (H + py. 

To proceed, we consider the special behaviour of the solution in the jump 
region. Let us first consider the case where the limit surface Sx given 
by Lemma [12] is not a vertical cylinder. Then it is a graph, which means 
that \Vui\ converges locally uniformly to something finite, and therefore 
that |Vnj| = e^Vi^l converges locally uniformly to 0. In the other case the 
surface Sx given by Lemma [12] is a vertical cylinder. We know from (|72p 
and (|75p that \VUj\ converges in I? to something finite. However this limit 
can only be zero since Ui — >■ U locally uniformly, and U is constant in the 
jump region (namely U = to on JCt by hypothesis). Furthermore, since the 
local uniform convergence of Ui —> to holds for the entire sequence i, we 
must have L 2 convergence of the entire sequence |VE7i|, which implies 

\VUi\ 0. 

□ 
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Proof of Proposition \22l Proposition [8] and Lemma [23] imply that each sur- 
face Sx in the jump region bounds a Caccioppoli set that minimises area 
plus bulk energy (g i3 — v % i>i\ Kij in A^ . To complete the proof of Proposi- 
tion El it remains to show that each surface £x in /Q is in fact a smooth 
MOTS. To proceed, we use the connection between parametric and non- 
parametric variational problems, that follows from the relationship between 
a function w G BV\ oc {Sl) and its subgraph 

(76) W = {(x, t) G n x R : t < w(x)}. 

In particular, let ipw denote the characteristic function of the subgraph (|76|) . 
Then Theorem 14.6 in [Gi] states 

(77) f ^l + \Vw\* = / \V(pw\. 

In (H5J) we established that at each point Xq E K, to , there exists a subse- 
quence ij and a function w G C ' a (Bn(xo)) such that 

Wij ->■ u> in C^Br^o)), 

where B R (x ) := Tn^ +1 (io), where 

(78) graph(^) = t Xo = exp" 1 (t Xo n B% xR (X ) 



Then Lemma [13] establishes that each surface T,x bounds a Caccioppoli 
set E in /Q that minimises area plus bulk energy (g lJ — ^D^Kij in /Ct , 
where by construction v is the outward unit normal vector to the relative 
boundary dE n ICt . 

Therefore, writing the Caccioppoli set E locally as the subgraph of w := 
exp~ 1 (i/)), we find from (|77p that it? minimises the functional 



J? r(x °\w):= Jl + \Vw\ 2 dx 

JB r (xo) V 



/■ pw(x) 

+ / trMKij{x,T) — Kij(x,T)i> l (x,T)i>i (x,T)d,Tdx 

Jb r (x ) Jo 

in Br(xo) := exp 9 (I?RXo))) whose Euler-Lagrange equation is the MOTS 
equation 



(79) div + (gV - PV) K 




and by construction v = — — The left hand side of (|79]) is an 



elliptic operator of the form 



2 



^ = a y (Vw) V.Vj-u; + Jl + | Vw\ 2 K 
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where 



o«(p) := f== [g 



1 ( _ p l p> 



y/i+w y 1 + \p\ 2 

Since w G C X ' a (B R {xo)), a ij G C°' a (B R (x )), Aw is strictly elliptic on 
Br(xo). Schauder theory then implies that w G C 2,a (B R (xo)), and by boot- 
strapping further we obtain w G C°°(B R (xo))- Using a suitable partition of 
unity, we obtain that each surface Sx is a smooth MOTS in /Cf . □ 

6. Existence of weak solutions 

In this section we use the normal vector field v of the hypersurface foliation 
of the jump region Kt Q from Proposition [8] to extend v = j^jj across the 
jump region, thereby constructing a globally defined normal vector field 
v. Existence of weak solutions is then proven by taking the limit of the 
translating graphs Sf, using Compactness Property 

Theorem 24 (Existence of weak solutions). Let (M n+1 ,g, K) be a complete, 
connected, asymptotically flat initial data set without boundary, that satisfies 
trjtfK > 0. Then for any nonempty, precompact, open set Eq C M with C 2 
boundary, there exists a weak solution of (**) with initial condition Eq. 

Proof. Let U be the limit of U e as given by (|32p . We construct the vertical 
cylinders := d{U < t} and T,f := d{U > t} with local uniform C 1,a 
bounds and unit normal vector field v with local C°' a bounds. Then using 
Theorem ([9]), we show that {U < t} minimises ^Tjj,p m (M\Eq) x R for each 
t, where v is extended in the jump regions tCt Q by the normal vector field v 
to the family of smooth MOTS {^x }xo£fct ' 

i) In the case where St = T,f, the surface St is constructed by fixing a 
point X = (xq, zq) € St and considering the sequence of times ti such that 
Xq G SJ. for each i. It then follows exactly as in the proof of Lemma [T2l that 
Sj. converges locally uniformly to St. Since Sf = S^" is a vertical cylinder, 
convergence holds for the full sequence, and the unit normal vector field v 



is equal to 



iwr 



ii) We use a slightly different pointwise approach to construct Sf and S^" 
when St 7^ S^". To this end, let Xq G S^ at a jump time tQ. Since there 
are only countably many such to, there exists a sequence of points X{ G S^ 
with t{ > to, satisfying lim X{ = Xq and lim ti = to. For i > 1 large 

enough, we can assume that S^ = Sj, and as above each surface piece 
S^ n Bft X ^(Xi) can therefore be written via the exponential map (denoted 
by the hat superscript) as the graph of a C 1,a function Wi over S^, where 

t ti :=e W -](t ti nB^(Xi)). 
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Now consider the sequence i>{ of normal vectors to at Xi. Since the 
C'i(Xi) are uniformly bounded in C 0,a , there exists a subsequence and a 
unit vector field v such that i>i . — > u uniformly on B n -^ 2 (Xi ) . Let T denote 

the hyperplane containing Xq and orthogonal to u(Xq). For i > 1 large 
enough, we can then write each surface locally as the graph of a C 1,a 
function vbi over T n B r ^ 2 (Xo). By Arzela-Ascoli, there exists a further 

subsequence and a C 1 '" function m:Tn -B^ +1 (Ai) such that 

wi^w inc 1 (rn^+ 1 (x i )), 

where Xo G graph u) and T = Txgraphtu. In order to recognise graph -u) as 
a piece of E^ and T as T^- Sj£, we consider a point Y G graphs. Then 

there exists a sequence Y G graph Wi C E^ such that Yj — )■ Y", and thus 
?7(Yj) = tj implies ?7(Y) = to, where U := U o exp . 

In order to obtain a contradiction, assume that Y G Ef Q . Then there 
exists 5 > such that B^ +2 (Y) £ E?^, however this contradicts the fact 
that Yj G T'ti for > to- Thus we deduce that graphs E E^ as required. 
The case where Xq G E ro for Ef. 7^ E^ follows analogously. 

In Proposition [8] we constructed a family of surfaces {E A - o }x <=jc t foliating 
the jump region ftt of U . This foliation has a C, normal vector field v, 
which extends the vector field of the surfaces Et and E^ as a calibration 
across the jump region at jump times to via the definition 



v{x) := < 



w 
\vu\ 


(x) 




if x 




at regular times t, 




V 






if x 


G /C ro 


at a jump time to, 




lim 

i— >oo 


w 

|W| 




if x 


G E to 


where Xj G E^ for Xj - 


■> X,tj /" t , 


lim 

2— >-oo 


VC/ 




if x 




where Xj G E^ for Xj - 


■> X,ti \ t . 


IWI 



This global interpretation of the normal vector field v in M\Eq x M means the 
functional Ju,v is well defined on M\Eq xi, and it follows from Compactness 
Property [9] that the sets {U < t} minimises Ju,v in M\Eq x K for each t. 
The result then follows from Lemma fl3l □ 

7. Applications of weak solutions 

In this section we highlight the natural applications of weak solutions of 
(**) to the existence theory for MOTS and to the theory of weak solutions 
of IMCF. 
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MOTS. The one-sided variational principal associated with outward op- 
timisation implies that the solution must jump at t = wherever the null 
mean curvature of So = 8Eq is strictly negative. Together with Proposition 
[22] this implies the following existence theorem for MOTS in initial data 
sets (M,g,K) containing an outer trapped surface So such that 6^ o < 0. 

Proposition 25 (Existence of smooth MOTS). Let (M n+1 ,g,K) be an 
asymptotically flat initial data set of dimension n + 1 < 7 satisfying trMK > 
0, and let Eq be any nonempty, precompact, smooth open set in M satisfying 
0q Eq < with respect to the unit normal pointing out of Eq. Then the level 
set d{u > 0} of the weak solution (U(x,z) = u(x),p) of (**) is a smooth 
MOTS in M\E . 

Proposition [25] highlights the natural application of the theory of weak 
solutions to (**) as a variational- type approach to constructing marginally 
outer trapped surfaces. Let k denote the largest eigenvalue of K with respect 
to g across M\Eq. Then in the special case where the outermost MOTS in 
M\Eq satisfies 

(80) \dE\ < \dF\ - nKC n+1 (E\F), 

for every closed MOTS dF in M\Eq, then the weak solution to (**) in 
Proposition [25] will jump to the outermost MOTS S = d{u > 0} at t = 0. 
In general, given any initial data, and any initial condition Eq satisfying 
0q Eq < 0, the surface can't jump beyond the outmost MOTS at t = 0, 
and the MOTS d{u > 0} is an inner barrier for the outermost MOTS S in 
M\E . 

We compare Proposition [25] to the following existence theorem combining 
[2] and [6], as stated in pp. 

Theorem 26 ([2] [6]). Let (M,g,K) be an initial data set of dimension 
n + 1 < 7 and let O C M be a connected bounded open subset with smooth 
embedded boundary dVt. Assume this boundary consists of two non-empty 
closed hypersurfaces d+£l and d-£l, possibly consisting of several compo- 
nents, such that 

(81) H d+n - tr a+n K > and H 9+n + tr 9+n K > 0, 

where the mean curvature scalar is computed as the tangential divergence of 
the unit normal vector field that is pointing out of Q . Then there exists a 
smooth closed embedded hypersurface S C O homologous to <9_f2 such that 
H-£ + tr^K = (where H% is computed with respect to the unit normal 
pointing towards d-£l). S is X-almost minimising in $7 for A = 2(n + 1)k, 
where k denotes the largest eigenvalue of K with respect to g across fL 

Here A- almost minimising in f2 means that the surface S arises as (a 
relative boundary of) a subset E C with perimeter S in Q such that 

(82) \dEnW\ < \dFHW\ + \£ n+1 (EAF), 
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for every FcSJ such that EAF CC W CC f2. A detailed analysis of such 
almost-minimising boundaries is carried out in [5]. We say that the set E is 
A-almost minimising on the outside/inside in Q if E satisfies (|82p for every 
F such that EAF CC W, where F C E, F D i£ respectively. 

Remark. Since minimising area plus bulk energy P is a stronger con- 
dition than A-almost minimising, weak solutions (U(x,z) = u(x),v) of (**) 
satisfy the following A-almost minimising properties, for A := nn, where k 
denotes the size of the largest eigenvalue of K on M\Eq: 

{%) The smooth MOTS d{u > 0} of Proposition[25lis A-almost minimising. 
(ii) The sets E± = {U < t} and {U < t} are A-almost minimising on the 

outside for each t > 0, t > respectively, 
(m) The surfaces £x foliating the interior tCt of the jump region are A- 
almost minimising in /Ct . 

IMCF. In this section we discuss the above results in the context of the 
work of Huisken and Ilmanen [10] on inverse mean curvature flow. In par- 
ticular, when applied to the special case K = 0, Definition [29] provides a 
new perspective on weak solutions to inverse mean curvature flow, and the 
work of sections [3] and [5] carries over to prove the analogous results for the 
jump region. 

The degenerate elliptic equation 

(*) div M = \VU\ 

describes inverse mean curvature flow of the level-sets of the scalar function 
u : M — > R wherever |Vu| 7^ 0. In [10] . Huisken and Ilmanen define a 
locally Lipschitz function u G C^ c (M) to be a weak solution of (*) with 
initial condition £^0 if Eq = {u < 0} and 

(83) J u (u) < J u {v), where J u {v) = j£{y) := / |V«|+«|Vu|, 

J A 

for every locally Lipschitz function v with {» / u) CC M\Eq, where the 
integral is performed over any compact set A D {u 7^ v}. They then showed 
that it follows that u is a weak solution if and only if the open set Et := 
{u < i] minimizes the parametric energy functional 

(84) J*(F) = \d*FnA\- [ |V«|, 

J FnA 

in M\E for each t > 0. 

Theorem 27 (Existence of weak solutions, [10J). Let M be a complete, 
connected Riemannian manifold without boundary. Suppose there exists a 
proper, locally Lipschitz, weak subsolution v of 183\) with a precompact initial 
condition. Then for any nonempty, precompact, smooth open set Eq in M, 
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there exists a proper, locally Lipschitz weak solution u of (*) with initial 
condition Eq, which is unique on M\Eq. 

In [9], Heidusch proved optimal Cj ' regularity for the level sets N t = 
d{u < t} and iV t + = d{u > t} of the weak solution. The theory of weak 
solutions to inverse mean curvature flow as laid out in |10j . however, does 
not include an analysis of the interior of jump regions. Applying Proposition 
[22] in this special case where K = 0, we obtain a foliation of the interior of 
the jump region {u = to} x 1 by area minimising hypersurfaces, a result 
which was left open in |10j . 

Corollary 28. Let u be the weak solution of inverse mean curvature flow 
given by Theorem \27\ At a jump time to, the interior JCt Q of the region 
{u = to} X R is foliated by smooth area minimising surfaces, each of which 
is either a vertical cylinder or a smooth graph over an open subset of /Q . 

We can then utilise the jump region hypersurfaces of Corollary [28] to 
present a new perspective on weak solutions of inverse mean curvature flow. 
In particular, by instead considering the weak solution to be a family of 
hypersurfaces one dimension higher in M x R, we can ask for the functional 
Jjj, defined by (|84|) . to be minimised everywhere in (M\Eo) x R, and obtain 
the following richer notion of weak solution. 

Definition 29 (Alternative weak formulation). Let u be the unique, locally 



Lipschitz weak solution to (■*•) on M\Eo given by Theorem 27, and define 
the locally Lipschitz function U(x, z) := u(x) on (M\Eo) x R. The weak 
solution to (*) is defined to be the pair (U,is), where v is a unit length, 

VU 

translation invariant extension of jy^j" * n ^ e 3 um P regions such that at 

each point x G JCt , v(x) is the normal vector to a C l,a hypersurface passing 
through x, which bounds a Caccioppoli set that minimises Jjj,u in K-t - 

We obtain the following weak existence result as a corollary of Theorem 



Corollary 30 (Existence of weak solutions). Let M be a complete, con- 
nected Riemannian n-manifold without boundary. Suppose there exists a 
proper, locally Lipschitz, weak subsolution of Ii83\) with a precompact initial 
condition. Then for any nonempty, precompact, smooth open set Eq in M , 
there exists a weak solution satisfying Definition{2Min M\Eq xR with initial 
condition Eq. 
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